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SUMMARY 


The application of the time-lag products correlation/ 
frequency analysis procedure to determine subcritical 
frequency and damping from structural response measurements 
made during flight flutter test of the B-l prototype air- 
plane is described. A detailed description of the analysis 
procedure is presented, and the test airplane and flight 
test procedures are described. Summary frequency and 
damping results are presented for six transonic flight 
conditions. Illustrative results obtained by applying 
various options and variations of the analysis method are 
presented for one flight condition. 


INTRODUCTION 


This study is one of several technology programs which 
use B-l flight test data and are being conducted in coopera- 
tion among the United States Air Force , ^Rockwell 
International, and NASA. Several features of the B-l flight 
flutter test program have enabled it to yield experimental 
data well suited for the evaluation of subcritical flutter 
response data analysis techniques. These features include 
extensive aircraft instrumentation, an onboard flutter 
excitation system, real-time telemetry monitoring, a data 
acquisition system, transonic/low-altitude flight capability, 
and a structural design based on flutter optimization 
analyses and flutter model tests. 

The purposes of this report are (1) to describe the 
important features of the flight flutter and data analysis 
programs, (2) to illustrate the application of the data 
analyses methods to B-l subcritical response data, and 
( 3 ) to show that the combination of the elements of the B-l 
flight flutter test program resulted in the generation of 
reliable subcritical frequency and damping information 
which allowed efficient and safe flight envelope expansion. 



Demonstration of the flutter stability of a high-performance aircraft 
to flight envelope limits is a critical development milestone. Because of 
the limits on accuracy of flutter analyses and flutter model tests, and the 
weight penalty involved if large margins of safety in design are used, flight 
flutter test demonstrations must be carried out on an actual prototype of all 
new high-performance aircraft. 

The flight flutter test must provide evidence of the existence of a 
flutter margin of safety above the maximum operating speeds, including over- 
speed increments. However, the risk of catastrophic failure is too great to 
allow flight near the predicted flutter point. Therefore, the flutter clearance 
of the flight envelope is accomplished by measuring the frequency and damping 
of the structural vibration response at subcritical speeds, starting with low 
speeds and gradually increasing speed up to a maximum speed which allows for 
the required margin on the predicted flutter boundary. At a given test condi- 
tion, these frequency and damping characteristics are extrapolated to the next 
higher speed or dynamic pressure condition to predict flutter stability for 
safe flight envelope expansion. 

This approach to demonstrating the flutter margin of safety makes the 
safety of the flutter clearance program heavily dependent on the determination 
of the response frequency and damping from the subcritical response measured. 
Consequently, the continued improvement of the methods used to determine 
frequency and damping from flight test response data is of great interest to 
all organizations connected with aircraft design and development. 

The basic intent of this study is to evaluate the basic elements required 
to obtain subcritical damping and frequency information. These basic elements 
include a method of exciting the structure in flight, instrumentation, test 
procedures, and analysis of the flight test subcritical response data. These 
evaluations will add to the technology base for future aircraft development. 
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entire time interval of a time history 
position transducer 
correlation lag number 
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maximum number of time lags 
number of cycles 

number of digitized data points or samples 

Newton’s metric force unit 

power 

dynamic pressure in kN/m^ 
quadrature spectra of response and force 
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autocorrelation function of response 

autocorrelation function of input force 
cross-correlation function of response and force 
cross-correlation function of force and response 


coherence function 

power spectral density of response 

power spectral density of force 

cross spectra of response and force 

cross spectra of force and response 

sample rate 

spectral 

sampling period in seconds per sample = 1/SR 
time 

time history of response 
time history of input force 

mean deviation = V' | V j - V [ where V is variable of 
j = 1 * 

interest, V is the average of n variables 
time length of correlation functions, seconds 
phase angle 

power spectral density, (units) ^ / RAP 

SEC 

frequency in radians/second 



FLIGHT TEST PROGRAM 


The B-l Flutter Test Aircraft 

General description. - The B-l aircraft, shown in figure 1, is a prototype, 
long-range supersonic banber with the capability of transonic and supersonic 
flight at high dynamic pressures. The aircraft utilizes a blended wing-body 
concept with variable- sweep wings, a single vertical stabilizer with a three- 
section (upper, intermediate, and lower) rudder, and horizontal stabilizers 

which provide both pitch and roll control. The variable-sweep (15 to 
67.5 degrees) wing, as shown in figure 2, is equipped with slats, spoilers 
(which also function as speed brakes) , and flaps which provide the aircraft 
with a highly versatile operating envelope. Canted vanes, mounted on each side 
of the forward fuselage, are part of the structural mode control system which 
reduces structural bending oscillations in the vertical and lateral axes. For 
flight flutter tests, a flutter excitation system (FLEX) was added to A/C-l. 

The aircraft is powered by four YF101-GE-100 dual-rotor augmented turbo- 
fan engines in the 30,000-pound-thrust class. The engines are mounted in twin 
nacelles below the wing and aft of the wing pivot points. For supersonic 
speeds, an air induction control system varies the internal geometry of the 
nacelle inlet ducts to maintain the required airflow to the engines for all 
flight conditions. 

Fuel is carried in integral tanks in the fuselage, wing carry -through, and 
wing outer panels. The fuel system is pressurized and inerted by nitrogen. 

Fuel transfer sequencing is automatic and provides center-of-gravity control. 

The aircraft has both in-flight and ground single-point refueling capabilities. 

Instrumentation . - Figure 3 shows the 43 accelerometers and strain gages 
located throughout the airframe to sense the response of the airframe and 
FLEX exciters. The outputs of all of these transducers are recorded in analog 
form on an onboard tape. The outputs of 21 selected transducers are telemetered 
for presentation on three strip charts at the ground station at Edwards Mr 
Force Base for real-time monitoring of flight safety and for evaluating the 
adequacy of frequency sweeps as they are completed. 

The primary response accelerometers used for data analysis are on the 
wing, horizontal tail, and vertical tail tips. Both left and right wings and 
horizontal stabilizers are instrumented to allow determination of the symmetry 
or antisymmetry of the response. The left wing, left horizontal tail, and 
vertical fin have response accelerometers near the leading and trailing edges 
of the tip chords to aid in identifying bending or torsion responses. Con- 
trol surfaces such as flaps, spoilers, and rudders are also instrumented. 

Other secondary components instrumented are the movable overwing fairing. 
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nacelle movable inlet ramp, structural mode control vanes, aft nacelle, and 
various fuselage locations. This array of instrumentation locations allows 
sufficient coverage of the vehicle to determine the subcritical response of 
the important components. 

The flutter excitation system. - The B-l flutter excitation system 
consists of five hydraulically driven, electronically controlled oscillating 
mass (or wand) inertial force exciters, mounted at the wing, horizontal tail, 
and vertical tail tips, as shown in figure 4, a. A photograph of an exciter 
mounted to a test stand is shown in figure 4, b. This inertial type of exciter 
offers certain advantageous characteristics, such as allowing the excitation 
forces to be independent of flight speed and wing sweep angle, and causing 
minimal aerodynamic interference on the surface. 

The FLEX control loop is shown in figure 5. Its major components include 
the copilot's control panel, the oscillator control, an overresponse safety 
trip system, the pilot T s emergency off switch, the hydraulic actuators and 
servo, a reference accelerometer mounted on the surface structure, and an 
accelerometer mounted on the wand. The wing and horizontal tail exciters are 
synchronized to operate either in phase for symmetrical excitation or out of 
phase for antisymmetric excitation. 

The copilot's FLEX control panel, shown in figure 6, contains switches 
and controls for selecting wing, horizontal tail, or vertical tail, setting 
inertial force levels, selecting symmetric or antisymmetric, setting the 
excitation frequency at a discrete value for dwells, performing automatic 
frequency sweeps, and manual system shutdown. In the automatic frequency 
sweep mode, the excitation frequency is swept at a logarithmic rate from 1 to 
65 Hz in approximately 55 seconds, as shown in figure 7. Reverse or down 
sweeps also can be made. The FLEX wand frequency is controlled by the oscil- 
lator signal. Due to the finite reaction time of the hydraulics to the sweep 
frequency command and a time lag introduced by the feedback rectifier, a 
frequency-dependent phase lag exists between the FLEX wand and the oscillator, 
as shown in figure 8. 

The response safety trip system continuously monitors the amplitude of 
13 accelerometers simultaneously. It provides automatic FLEX shutdown in the 
event that any one accelerometer indicates that wing, horizontal tail, or 
vertical tail is responding in excess of predetermined response limits. Auto- 
matic FLEX shutdown will also occur if the shaker forces should exceed pre- 
determined limits. The cutoff response and force levels are chosen to be 
well below the amplitudes that could cause any structural or FLEX damage. In 
the event of failure of the hydraulic supply, the exciter wand is hydraulically 
locked to prevent the wand from responding to airframe motion. 

The hydraulic actuator consists of two pistons operating in parallel, to 
directly drive the wand. Movement of the wand also drives a position indicator 
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(LVDT) which controls operation of the mass assembly about zero. An electrically 
operated hydraulic servo valve ports oil pressure to either side of the piston 
when commanded to move plus or minus by the electronics controller as a func- 
tion of frequency. 

Oscillation of the wand relative to the structure is sensed by the dif- 
ference in output from two accelerometers. One accelerometer is on the wand, 
and the other is a reference accelerometer on the aircraft structure, imme- 
diately adjacent to the wand. The wand accelerometer and the reference 
accelerometer are at the same chordwise location on the aircraft surface. The 
chordwise location selected is the unique location at which the wand excitation 
force acts on the structure. At this location there is an inertia force, but 
no moment, due to relative motion between the wand and the aircraft structure. 

The inertia force and moment of the oscillating wand are reacted at the wand 
pivot point by a force and a moment. However, this combination of force and 
moment can be considered to be replaced by an inertia force alone acting at a 
location which gives the same moment about the wand pivot point. The location 
of this inertia force is calculated from the force and moment reactions at the 
wand pivot. Figure 9, a, represents the wand with a small rotation, $, relative 
to the aircraft structure. The vertical force at the wand center of gravity is: 

F Cg ’ n ^g 

where 

y = distance from wand pivot to wand c.g. 
eg ^ & 

m = mass of wand 

The vertical reaction at the wand pivot is: 

F — p 

Pivot eg 


The moment about the wand pivot is: 


Mpivot 

where 

I = wand inertia about pivot 
p 
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The location at which the reaction force and moment on the structure can 
be replaced by a vertical force alone is determined by dividing the wand 
inertial moment by the inertial force to yield 


X 


e 



where 


X 0 = distance from pivot to effective force location 

An effective mass m at x can be defined such that its inertia force will 

0 0 

equal the wand inertia force and its moment about the wand pivot will equal 
the inertia moment of the wand. Equating the moments about the pivot due to 
the wand and due to the effective mass results in 


2 " 

m x 3 

e A e 


- V 


so that 



Substituting the preceding expression for x and S = my , the effective 


mass given is 


eg 



The wand inertia force then is given by 


F 

FLEX 


m Aa 
e 
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where Aa is the difference between the wand accelerometer and the reference 
accelerometer output. Any overall motion of the aircraft does not affect Aa 
since the wand and reference accelerometers are very close together. Any 
relative motion between the wand and the structure causes a nonzero Aa. 

Lightweight wands weighing 137 newtons (31 pounds) with an equivalent 
point mass of 79.5 newton (18 pounds) were used for early flight flutter tests. 
The response obtained with the light wands sometimes was not large enough to 
provide consistent danping and frequency data when the aircraft was subjected 
to turbulence and buffet. Heavier wands were installed for later flight 
flutter tests to obtain larger responses relative to the responses to tur- 
bulence and buffet. The heavier wands weigh 184 newtons (41.3 pounds) and 
have an equivalent point mass of 123.8 newtons (27.8 pounds). Figure 10 shows 
the input force as a function of frequency for the wings, horizontal tails, and 
vertical tail for the lightweight and heavyweight wands. Ideally, the force 
should be constant at all frequencies; however, figure 10 shows that the force 
amplitude was a function of frequency. The FLEX force increased with fre- 
quency from zero to about 10 Hz due to displacement limitations of the wand 
travel. From about 10 to about 35 Hz, the force level was fairly constant. 
Above 35 Hz, the force decreased due to the increasing phase lag of the wand 
relative to the oscillator command. This latter rolloff is more pronounced 
with the heavyweight wands. The vertical tail wand had the lightweight version 
for all flights. The higher force level for the vertical tail on flight 1-50 
compared to flight 1-13 is due to a higher amplitude setting on the copilot's 
FLEX control box. 


Flight Tests 

Test conditions . - Figure 11 shows the six flight-test conditions selected 
for analysis in this report. The test conditions consist of three altitudes 
at mach 0.90 and three altitudes at mach 0.95. The test points are arranged 
so that subcritical frequency and damping trends at constant mach number, 
constant dynamic pressure, and constant altitude can be determined by appro- 
priate plots of the test data. Test dynamic pressure increments between points 
at constant mach number were chosen to be small enough so that decreases in 
damping with increasing dynamic pressure should be safely detectable. FLEX 
frequency sweeps were performed, and 5 minutes of random excitation response 
data were recorded at each flight condition. The maximum buffet response of 
the flight envelope was encountered at these two mach numbers and increased 
in amplitude with decreasing altitude. 
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Test procedures. - The flutter tests were performed with, the wings at 
a 65-degree sweep angle. The aircraft was trimmed at the proper mach number 
and altitude for the chosen test condition. Five frequency sweeps were per- 
formed for each test configuration. All data presented in this report are for 
the stability and control augmentation system and the structural mode control 
system deactivated. The frequency sweeps performed were: 

(1) Wing symmetric 

(2) Wing antisymmetric 

(3) Horizontal tail symmetric 

(4) Horizontal tail antisymmetric 

(5) Vertical tail 

The sweeps covered a frequency range of 1 to 65 Hz and were about 
55 seconds in duration. The resulting data were analyzed in three different 
time frames, as shown in figure 12. Real-time monitoring of the response and 
force data involved observation of analog strip charts in the Mission Control 
Center at Edwards Air Force Base. The strip chart monitor determines the 
successful completion of a FLEX sweep by observing such characteristics as 
adequate signal- to-noise ratios, response amplitudes not causing automatic 
FLEX shutdown, visual judgment of adequate damping in all modes, atmospheric 
turbulence not being excessive, and proper functioning of the accelerometers. 

If the monitor observes any potentially dangerous response, he can immediately 
inform the pilot by direct radio communication to deccelerate. Near real-time 
analysis uses a minicomputer in the Mission Control Center to digitize the 
data as a sweep is performed for subsequent data analysis. The data are 
analyzed using a time-lag products correlation/frequency analysis method 
described later. The frequency and damping of the modes contained in one 
chosen frequency range are calculated from the gain or transfer function, 
usually within 20 minutes after the completion of the sweep. Postflight data 
analysis either uses the minicomputer to complete the analysis of all modes 
analyzed during the flight, or employs the onboard tape recording of the sweep 
data for offsite time-lag products correlation/frequency analysis using a large 
storage IBM computer to calculate modal frequency and damping. 
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SUBCRITICAL RESPONSE ANALYSIS 


Time-lag Products Correlation/Frequency Analysis Procedure 

The primary procedure used to obtain the subcritical response from B-l 
flight test data involves the time-lag products correlation/frequency analysis 
method. This method derives the time domain autocorrelation and cross- 
correlation functions for one pair of time history records, and estimates 
spectral density and transfer functions from these calculations using digitized 
data. Modal frequency and damping values are calculated from the various 
spectral density and correlation functions. Coherence estimates are also 
calculated. The procedure also includes an optional smoothing technique to 
improve the spectral estimates by the use of a "Hanning window," and an 
optional technique to calculate modal frequency and damping from the transfer 
function which uses an iterative, optimal, least squares complex curve fit. 

The analysis procedure first plots the digitized data in the form of time 
histories of the chosen response channel and the corresponding FLEX force 
channel, or the oscillator signal (which controls the FLEX frequency). The 
time histories are plotted without and then with the mean amplitudes removed. 

The procedure then generates five types of data representations from which 
the modal frequency and/or dampings may be manually calculated. 

The first type of data representation is the correlation functions. The 
autocorrelation function tends to identify a sinusoidal or periodic function 
buried in a single noisy time history, while the cross-correlation function 
identifies periodic functions which have common frequencies in two different 
time histories. Autocorrelation and cross-correlation functions are obtained 
by multiplying two time histories point by point and summing these products. 
Then, one of the waveforms is shifted or lagged by some time increment, and 
the sum of the point-by-point multiplication of the two wave forms is again 
calculated. For responses with stable damping, the greater the second time his- 
tory is time shifted, the smaller the sum. The correlation function of the two 
time histories is obtained by plotting the sums of the time-lag products against 
the time shift or lag number. The autocorrelation multiplies a single time his- 
tory and its identical time history, while the cross-correlation multiplies two 
different time histories. The autocorrelation and cross-correlation functions 
are calculated by the classical Backnan and Tukey (ref. 1) time-lag products 
correlation method. If the input (force) time series is called y (t) , the cor- 
responding autocorrelation function is denoted by Ryy (mT) and is given by 

N-m 

R (mT) = r^— ^ ' Yv ' Yv , m = 0, 1, 2...M 
yy v N-m 7 K 7 K+m’ 

K=1 


11 


where 


N = total number of digitized data points 
m = correlation "lag number" 

M = maximum number of time lags (= 500 for this report) 

T = sampling period in seconds = 1/ (sample rate) 

This calculation is performed for positive lag numbers only because Ryy is an 
even function ; i . e . , 


R (mT) = R (-mT) 

yy yy 


The output time series (response of the structure) is denoted by x(t). 
Its corresponding autocorrelation is denoted by R (mT) and is given by 

XX 

N-m 

R (mT) = / x T/r • x., , m=0, 1, 2,...M 

xx J N-m K K+m’ 

K=1 

Rjo^- is also an even function. R^ and Ryy are plotted as a function of lag 
number. This may be converted to time by knowing the time shift per lag 


At 

lag 


_ 1 _ 

SR 


T 


where SR is the sample rate of digitization in samples per second and At is 
the time shift of each lag. For a sample rate of 360 samples per second, each 
lag would be 


4t * If ’ 3So seConds 


or 360 lags = 1 second. 
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The cross-correlation functions are then calculated and plotted in the 
same form as the autocorrelation functions. The cross-correlation function is 
denoted R-^Cmt) and is given by 


R (mT) 
xy^ J 


N-m 


N-m *K * y K+m’ m °» 1 » 2, “ ,M 


K=1 


and the cross -correlation function R 
by 


yx 


(mT) between input y(t) and output x(t) 


N-m 


R (mT) 
yx 


■(S £ Y K ■ =W 1. M 


K=1 


These functions are neither odd nor even but satisfy the relation 


R (mT) = R (-mT) 
xy yx£ 


This property is used for calculating the cross-spectral functions later. 

The plots of the autocorrelation function of the response and the cross - 
correlation functions may be used directly to estimate the modal frequency and 
structural damping, g, for a single-degree-of-freedom system by treating the 
functions as a free-vibration decay with an initial displacement, as mentioned 
in ref. 2. Thus, the frequency is determined by counting the number of cycles 
occurring in 1 second, and the structural damping is calculated as the log 
decrement (figure 13, a) or 



n 


for g £ 0.30 


where 


n = number of cycles 


A q = amplitude of initial peak 

th 

= amplitude on the n peak 
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However, for analysis of a time history containing more than one mode, deriving 
the frequency and damping of a single mode would be very difficult, if not 
impractical , because the correlation function would contain information from 
all of the modes. Elimination of the unwanted mode may be accomplished for a 
swept -frequency, forced response by time windowing; i.e., by simply analyzing 
only that portion of the time history containing the mode of interest. How- 
ever, this has limitations due to minimum time requirements for accurate data 
analysis. A problem with using the autocorrelation of the response for a direct 
frequency and damping determination is that any noise response in the time 
history (such as caused by turbulence or aircraft buffet) that contains fre- 
quencies in the realm of the mode frequency of interest will distort the 
function. However, this effect is minimized when using the cross-correlation 
function because the process of cross -correlation emphasizes response data 
which have a common frequency as the force. In general, correlation functions 
were not usually used to directly determine modal characteristics due to the 
reasons already mentioned. However, they are used to generate power spectral 
(frequency domain) density and transfer functions where the modal frequency 
and dampings can be determined. 

The next step in the analysis procedure is to calculate the numerical or 
finite, discrete Fourier transform of the autocorrelation functions to derive 
the power spectral densities. The Fourier transform of the response is denoted 
as Sxxfjuj^) = Re[F(Rxx(mT) )] • Because the correlation functions are of finite- 
time duration, truncation errors in the Fourier transform may be introduced if 
the amplitudes of the functions at the end points are not small. This effect 
may be partially compensated for by introducing a two -point correction (ref. 3) 
for the end points with equal weights of one-half. Applying the end-point 
corrections, the Fourier transform of Rxx(mT) becomes 




= T 


R (0) 

xx 


M-l 


+ 2 


R (mT) 
xx 


cos 


mm, T 
n 


R (mT) 
xx 


cos 


NL^T 


m=l 


The power spectral density of the response is then plotted from which fre- 
quencies and dampings of modes may be manually calculated by the half-power 
point method shown in figure 13, b. This type of data display does not take 
into account any variation in force amplitude as the sxvept frequency changes 
nor does it minimize response caused by turbulence or buffet noise inputs. 
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Likewise, the autospectral density of the input force is 



R (0) 
yy 


M-l 

2 Xrf F yy C mT ) cos mo^T + R^CMT) cos Mu^T 
m=l 


The cospectral and quadrature spectral functions are also used to cal- 
culate response frequencies and dampings. The cospectra is the component of 
the total response vector which is in-phase with the force vector, and the 
quadrature spectra is the out-of -phase component. The modal peak frequency is 
most accurately determined from the cospectra, and the damping calculation is 
made using the quadrature spectra, as shown in figure 13, c. 



where 


The cospectra, denoted as C , and the quadrature spectra, denoted as 
are given as 

yay = y F{ p c (mT:i}] 

yoy = RjF{p q CmT)}], h = 1, 2. . . ,M 
p CmT) = 1/2 [R CmT) + R (ml)] 

^ a y y a 

p q (mT) = 1/2 [R^ (mT) - R^mT)] 


The Fourier transform with end-point correction factors gives 


V j V 


M-l 


V 3 "h> ■ 


p c C0) + 2 ^ P c (mT) cos moj^T + P c (MT) cos Mi^T 

(an even function) 
P q (mT) sin moj^T + p q (MT) sin Maj^T 
(an odd function) 


m=l 


M-l 


m=l 
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The coherence function of the response and the force is also calculated. 
Coherence is a measure of the degree to which two time series x(t) and y(t) 
are related by a linear model, as shown in figure 14. In this case, y(t) is 
the FLEX force input, X*(t) is the airframe response from the FLEX only, and 
N(t) is extraneous noise inputs such as turbulence and buffet which are 
uncorrelated with y(t) . x(t) is the total response of the airframe, and G(jw) 
is the linear model or the transfer function. Coherence values range between 
zero and 1.0. A high coherence (say, greater than 0.5) means that the response 
and the force at a given frequency are related. The coherence function, 
denoted by (oj) , is given by 

2 C w H “ ) * < 4«“ ) 

R (o>) = 7 T 

S (]Oj)S (jco) 
xx yy 


and is plotted against an The function is used to help interpret the validity 
of a spectral peak; i.e., if a peak frequency has a corresponding high 
coherence, then it is assumed that the mode was predominantly excited by the 
FLEX force rather than by buffet or turbulence noise; therefore, the peak has 
minimum distortion and should yield valid damping and frequency values. 

The frequency spectra of the noise may be obtained by knowing the spectra 
of the total response and the coherence function for the frequency range of 
interest. The noise spectra, denoted as S^Cjo)) , is given by 

S (jto) = S (jaj)[l - R 2 (oj) 1 
nn J xx 

The gain or transfer function is then calculated. This may be regarded as the 
frequency response function relating the input y(t) to the output x(t) . The 
gain is denoted by |G(jw)| and is given as 


I G(jw) | 



The gain frequency spectra are then plotted 
and damping are manually calculated ; as shown in 


from which the modal frequency 
figure 13, d. 
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The phase angle of the transfer function (or the phase between the FLEX 
force and response) is calculated and plotted. The phase is denoted by 0(jco) 
and is given by 


0(jaj) = tan 

\ xy 

The modal damping is calculated by using the slope of the phase curve, as is 
shown in figure 13, e. 

The Hanning filter or window is an option which may be applied to smooth 
spectral estimates. This filter is described in ref. 1 and 4. This program 
uses weighting coefficients of one-fourth, one-half, and one-fourth as the 
Hanning filter weights and the coefficients are applied for S^, Cxy, and Qxy> 
respectively. 

Another option that may be applied is a least squares curve fit which 
tries to match the gain and phase frequency responses. This method fits a 
mathematical transfer function to autospectral and cross -spectral density 
estimates of data by an iterative process. The analyst must specify the order 
of the numerator and denominator of the mathematical transfer function. From 
the mathematical transfer function fit, the modal frequency and damping are 
derived. 


Subcritical Response Analysis Variables 

Accurate data analysis requires careful selection of the numerical values 
for the analysis variables. Improper combinations of sample rate, the time 
length of the data sample, number of time lags, and frequency resolution can 
yield frequency and damping values with errors. 

Sample rate. - The minimum sample rate must be at least twice the highest 
frequency of interest to prevent aliasing. (A sample rate of five to six times 
the highest frequency is usually desirable.) The maximum sample rate which 
was used for the data analysis in this report was 360 samples per second with 
a low pass filter of 62.5 Hz. The highest frequency of interest was 60 Hz. 

This allowed a maximum-sample-rate-to-maximum-frequency ratio of 6. However, 
lower sample rates for analysis of low-frequency, lowly damped modes were 
required to prevent truncation of the correlation function, as discussed later. 
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This is accomplished by decimating the original data; i.e., using only every 
other digital point of the original 360 samples per second or an integer 
divisor thereof. Sample rates as low as 90 samples per second were used for 
the low-frequency ranges. 

Time length of data sample . - For swept sine response data, the time 
length of the data sample required for analysis is largely determined by how 
much of the time history contains the frequency range of interest. However, 
the time length of the data sample affects the accuracy of the correlation and 
spectral density calculations. Truncation of the swept sine time history must 
be at the points where the frequencies of interest are not present and which 
allow sufficient time history duration; otherwise, truncation errors may be 
introduced into the correlation and spectral functions. 

The total number of data points or samples N is the product of the time 
length of the data sample L times the sample rate (SR) or 

N = L(SR) = L/T 


To allow a minimum of two degrees of freedom to exist for a chi-squared 
distribution of the spectral estimates, the ratio of the number of data points 
(N) and the total number of time lags (M) should be 


N 

— > 
M “ 


2 


A ratio of approximately 20 is usually preferred. The data analysis in this 
report used sample lengths of from 12.5 to 25 seconds, yielding a maximum of 
18 chi-squared degrees of freedom for 500 lags. 

Number of time lags. - Five hundred time lags were used for data analysis 
in this report. The number of time lags and the sample rate determine the time 
length of the correlation functions. If a given frequency contained in a 
response or cross -correlation function has not fully decayed at the last lag, 
truncation errors are introduced into the frequency spectra, and the peak 
frequency and damping will contain errors. The end-point corrections previously 
discussed partially correct for this. This truncation may further be minimized 
by increasing the time length of the correlation functions by increasing the 
number of time lags and/or decreasing the sample rate. The time length of the 
correlation functions is given by 


(sample rate) 
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Suitable M and T variations are made until the decay amplitude of the 
desired frequency contained in the correlation function is small enough at t 
so that the effects of truncation are negligible. Care must be taken that the 
sample rate does not become so slow that aliasing could occur, or that the M 
becomes too high compared to N. 

An estimate can be made of the correlation time required for the function 
to decay to a small amplitude at the truncation time. To provide initial 
choices of M and T, x is estimated by 


t = MT = 


1 

fgTT 


In 


A 

o 

A 


T 


where 

f = estimated frequency 
g = estimated damping 

A q = maximum correlation function amplitude 
A = correlation function amplitude at t 

T 

For example, for an estimated f = 10 Hz, g * 0.06, and A /A = the minimum 
correlation time would be 


T 


1 

(10) (0.06)tt 


In 10 


= 1.126 seconds 


For a chosen sample rate of 360 samples/second, the required minimum number of 
lags would be 


M>ifS 

M * dm i 458 lags 
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If the analysis results prove to be inadequate, a higher M and/or a lower 
sample rate may be chosen and the analysis repeated. 

Corrections for damping errors introduced by truncation are discussed in 
ref. 5 for a boxcar window, and in ref. 4 for a Hanning window. 

Frequency resolution. - The incremental frequency chosen for calculating 
the spectral functions was 

DF > ^ (Hz), N < 10,000 


For NT = 25 seconds, DF = 0.04 Hz. To minimize any spectral errors caused by 
spectral window smearing, the desired minimum number of lags for a given DF is 


M = 


1 

2 DF T 


M < 500 


This condition cannot always be met due to storage limitations of the computer 
for a chosen DF and T. The incremental frequency for the data analysis in this 
report varied from 0.04 to 0.08 Hz. 


SUBCRITICAL FREQUENCY AND DAMPING RESULTS 
Summary Results 

The variations of the predicted subcritical frequencies and 
dampings with dynamic pressure for the six flight conditions 
listed in figure 11 are presented in figures 15 and 16 . These 
results were obtained using FLEX excitation and are values 
averaged from several calculations using different options and 
variations of the time-lag products procedure. 

Illustrative Analysis Results and Interpretation Techniques 

Results from FLEX and random excitation (buffet and turb- 
ulence) obtained by using different options and variations of 
the analysis procedure are presented in tables I through V, 
and in figures 17 through 24 for a Mach number of 0.95 and an 
altitude of 1,79^ meters (5,890 feet). This flight condition 
was chosen because the highest level of airplane buffet and 
turbulence response occurred here, thus providing response data 
that could be used to test the ability of the analysis method 
to separate modal responses excited by the FLEX system from 
responses induced by random excitation, and also to evaluate 
the capability of the analysis method to derive accurate modal 
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response from purely random excitation. Five modal responses were studied- 
wing symmetric second bending, the horizontal tail symmetric second bending, 
and three vertical tail/empennage modes. These modes were identified by 
comparing the frequencies of the peak responses with the flutter analysis 
results. Each mode presents a different type of analysis problem, as shown 
in figures 17 and 18. Figure 17a, shows the transfer functions (gain) 
generated by the FLEX frequency sweeps, and figure 17, b, shows the power 
spectra of the response from random excitation only. The wing response is an 
example of a generally "clean," well-separated mode with a moderately high 
FLEX response to random response ratio. The horizontal tail response is an 
example of a noisy mode with a relatively low-level of response amplitude. 

The vertical tail/empennage responses are examples of three modes which are 
closely spaced together in frequency and are affected by buffet forces. These 
choices were also influenced by the fact that the complete aircraft flutter 
analysis predicted that the wing, horizontal tail, and two of the vertical 
tail /empennage modes are important to the flutter characteristics of the air- 
plane. The portions of the time histories containing the frequency ranges 
of interest of the forced response, FLEX force, oscillator signal, and response 
to random excitation are shown for each mode in figure 18, a through c. 

Figure 19, a through g, shows the data analysis results for the wing mode 
excited by the FLEX system, which is considered a clean, well -separated 
response. Figure 19, a, shows the cross -correlation function decay calculated 
for 500 lags using a 25-second time history length sampled at 360 samples per 
second. The ideal function should appear as a log decrement decay with a 
constant frequency; however, the cross -correlation function is at best a rough 
decay with more than one frequency. This is due to the fact that a 25 -second 
time segment of the frequency sweep was used which contains three other wing 
modes which are all present in the correlation function. Therefore, this 
correlation function does not lend itself for a direct calculation of the log 
decrement damping without further windowing. 

Figure 19, b, shows the power spectra of the response (PSD) derived from 
the autocorrelation function for a frequency range of 36 to 66 radians per 
second (5.7 to 10.5 Hz). Since this calculation does not account for the 
relationship between the force and the response, it is most affected by buffet 
and turbulence noise inputs. A comparison of this function and the transfer 
function of figure 19, f, shows that the two minor peaks at 47.5 and 
55.75 radians per second on the sides of the peak at 51.5 radians per second 
are more pronounced on the PSD than the transfer function, which indicates 
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they are largely caused by noise inputs. This means that the PSD peak 
generated from the autocorrelation function is distorted more by noise than 
is the transfer function peak calculated from the cross -correlation function. 
This is demonstrated by the fact that the PSD required a much greater extrapola- 
tion of the sides of the peak (figure 19, b) to the half-power point than the 
cross-correlation peak required. (Extrapolation of the sides of the peak as 
shown in the figure was a standard procedure when the peak was judged to con- 
tain distortions). However, the half-power point PSD damping was g = 0.126 
and the peak frequency was 8.08 Hz, which is still in good agreement with the 
results of the other data displays which follow. 

Figure 19, c, shows the same wing data displayed in the form of the 
cospectra (the in-phase component of the total response) which gives the most 
accurate determination of the peak frequency, which was 8.2 Hz. 

Figure 19, d, displays the same wing data as the quadrature (out-of-phase 
component of total response). The function is a smooth curve; however, the 
shoulder peak to the right of the peak frequency of 51.5 radians per second is 
not well-defined due to the two places where the slope of the curve changes. 
Choosing the first shoulder on the right and the obvious shoulder on the left, 
the damping is g = 0.135. 

Figure 19, e, shows the coherence spectra (coherence is a measure of the 
degree to which the response and force time series are linearly related) . A 
large coherence value at a given frequency (say, greater than 0.5) may be 
interpreted as meaning that most of the calculated response spectra were 
directly caused by the force and only slightly affected by extraneous noise 
such as buffet or turbulence noise. Figure 19, e, shows that the coherence 
in the frequency range of the half-power point of the wing mode varies from 
0.6 to 0.9, which is interpreted as meaning the FLEX response is well -separated 
from any response caused by noise. Therefore, the modal data should be valid. 

Figure 19, f, shows the transfer function spectra for the wing mode. The 
figure shows a clean, symmetric, well-defined transfer function (gain) with the 
peak frequency of 8.16 Hz and half-power damping value of 0.119. 

Figure 19, g, displays the same wing data as the phase angle between the 
response and the FLEX force. The upper figure shows a reasonably smooth curve 
with a discernible slope for the left frequency range of the mode. The dashed 
line is an estimate of the slope and yielded a damping of g = 0.146. 
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Averaging the frequency and damping values from all displays for each 
basic analysis method yields an average frequency of 8.14 Hz and an average 
damping of 0.132, with corresponding mean deviations of 0.043 Hz and 0.009 for 
the damping. These deviations are considered highly acceptable. 

Figure 20, a through g, shows the data results for the horizontal tail 
mode, which is considered a noisy mode with a low response amplitude. The 
correlation function (figure 20, a) used 500 lags on a 25-second time history 
length sampled at 360 samples per second. 

These analysis results offer interpretation problems. For example, the 
half-power point amplitude on the power spectra of the response (figure 20, b) 
and on the transfer function (figure 20, f) is below three peaks. The largest 
peak is obviously the mode's peak frequency, but it is not obvious whether the 
sides of the major "clean" peak should be extrapolated downward to the half- 
power point (dashed lines) to calculate the damping, or whether the entire 
width of all three peaks together should be considered to belong to same mode. 
The most conservative approach would be to use the lowest damping value of the 
extrapolated peak, yielding g = 0.06. Using the entire width covered by all 
three peaks gives g = 0.152. One aid in deciding which value to use is to 
compare them with values obtained from other displays of the same data. How- 
ever, the other displays also offer the possibility of two different damping 
values. Averaging both sets of values from all four displays gives the extrapo- 
lation technique g * 0.0716 and, using the width of all three peaks, g = 0.132, 
with a common average frequency of 16.86 Hz. The corresponding mean deviations 
for the average damping values are 0.0072 and 0.027, respectively, which favors 
the conservative extrapolation technique. 

Figure 21, a through g, shows the analysis results for three closely 
spaced vertical tail/ empennage responses which are separated by only about 1 Hz. 
The transfer function of figure 21, f, adequately separated the three responses 
and yielded well-behaved functions. The other data displays, however, do not 
always clearly separate the first two peak responses and inadequately define 
the third peak, which had a low response amplitude. For example, the power 
spectral density of the response of figure 21, b, shows the two major peaks 
are poorly separated and therefore mutually distort their spectra shapes between 
them, as evidenced by their lack of spectral symmetry about their respective 
peak frequencies. As a standard procedure in cases such as this, the dampings 
are calculated by using only the half of each peak which is relatively 
unaffected by the other mode, and then doubling the frequency band at the 
half-power point by assuming the other side of the spectra is a mirror image. 
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Thus, the left side of the peak at 10.62 Hz (66.8 radians per second) was 
used and doubled for the first peak. The right side of the peak at 11.5 Hz 
(72.25 radians per second) was used and doubled for the second peak. The 
third peak at 12.73 Hz (80 radians per second) is not defined at all on fig- 
ure 21, b. The quadrature spectra (figure 21, d) adequately defines the 
response for the first two modes but not the third, in spite of the fact that 
the coherence (figure 21, e) shows high values for the third peak. 

Averaging the different displays for figure 21, a through g, yields for 
the first mode a frequency of 10.67 Hz and g = 0.058, with a mean deviation of 
0.033 Hz for frequency and 0.007 for damping. The analysis for the second 
mode yields average values of 11.64 Hz for the frequency and g = 0.084, with 
corresponding mean deviations of 0.097 Hz and 0.028. The high mean deviation 
of the damping for this mode may warrant only the use of the transfer function. 
The accuracy of the third mode values are in question because only the transfer 
function of figure 21, a, adequately described the mode. The test data yielded 
values of frequency = 13.13 Hz and g = 0.068. 

These results are summarized in tables I through V. They indicate that 
for a variety of modes the time-lag products correlation/ frequency analysis 
procedure provided reasonable correlation and spectral functions when excited 
by the FLEX system. These functions provided reliable frequency and damping 
results. 

Figures 22, 23, and 24 show the response of the wing, horizontal, and 
vertical tails to random excitation for a 25-second time history. The flight 
condition is the same as the previous comparisons which used the FLEX force 
excitation. 

Figure 22, a, shows the autocorrelation functions in the time domain 
(360 lags = 1 second) for the left and right wingtips. The dominating fre- 
quencies are 3 Hz (wing first bending), 8 Hz (wing second bending, visible only 
for the right wing), and 27.5 Hz (wing second torsion). Figure 22, b, shows 
the power spectra of the autocorrelations of the left and right wingtips along 
with the phase angle between the two. It was expected that the phase angle 
plot would indicate if a particular peak was a symmetric or an antisymmetric 
mode. Two distinct modes appear on both wing spectra. The first peak at 
6.6 Hz appears to be too low in frequency to be a wing mode, but it could be 
a fuselage mode. The second peak at 8.1 Hz does correlate in frequency and 
damping with the wing symmetric second bending mode obtained from the FLEX 
excited response; however, the phase plot shows a 160-degree phase angle, 
indicating an antisymmetric mode. However, the antisymmetric mode was measured 
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at a higher frequency of 8.99 Hz (55.5 radians per second) from the FLEX 
excited response and appears as a shoulder at 55.5 radians per second in 
figure 22, b. Thus, it appears that a reliable determination of symmetry 
from this method is impractical. Figure 22, c, displays the spectra with a 
narrower frequency window which eliminates the peak at 6.6 Hz. This yields a 
more coarse frequency resolution. However, the shoulder at 55.5 radians per 
second (8.83 Hz) now appears as a distinct peak and does correspond in fre- 
quency to the antisymmetric mode. Windowing caused the damping of the mode 
at 8.1 Hz to decrease by 50 percent compared to the nonwindowed mode. 

Figure 23, a through c, shows the horizontal tail response from random 
excitation. The response levels of the second bending mode are so small that 
the spectra of figure 23, b and c, appear as noise, and no useful information 
can be obtained from them. 

Figure 24 shows the vertical tail response from random excitation. The 
two distinct modes yielded frequencies and damping values (figure 24, b and c) 
in rough agreement with the FLEX excited results using timelag products auto- 
correlation. These data results using random excitation appear in tables I 
through V. 

Figure 24 shows the autocorrelation of random response data did give some 
reasonable answers for peak responses with relatively low damping. However, 
random excitation did not adequately excite all of the important modes to 
sufficient amplitudes for accurate determination of the peak frequency and 
damping using this analysis procedure on the 25-second time sample. Definite 
determination of symmetric or antisymmetric response was also found to be 
impractical. The random excitation results may have improved if longer time 
data samples were analyzed or if the vehicle could have encountered higher 
turbulence amplitudes. 

In the appendix are other variations of the analysis procedure. Included 
are using a Hanning filter, using the oscillator signal instead of the FLEX for 
cross-correlation, employing a curve fit to estimate frequency and damping, 
and time and frequency windowing. These variations did not yield significant 
improvements for the cases studied. 
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CONCLUDING REMARKS 


The time-lag products correlation/frequency analysis procedure is des- 
cribed, and results from application of this procedure, including many of its 
variations and options, to B-l prototype airplane flight flutter test data are 
presented. Subcritical frequency and damping trends obtained from analysis of 
response data, excited by onboard inertial shakers using sinusoidal frequency 
sweeps, are presented for six transonic flight conditions where buffet and tur- 
bulence effects are appreciable. Examples of the analysis procedure applied 
to peak responses presenting various analysis problems are discussed. The 
examples include modes with noisy peak responses and low response amplitudes, 
closely spaced modes affected by buffet and turbulence inputs, and clean, well- 
separated modes. The results showed that the analysis procedure adequately 
described the peak responses excited by the onboard shakers. This enabled 
reliable frequency and damping values to be obtained. The analysis procedure 
was also applied to random excitation responses. The results showed that 
reasonable frequency and damping information could be obtained by random 
excitation for some of the modes which were adequately excited and had a 
relatively low damping. However, in some cases reliable frequencies and 
damping of important modes could not be obtained from response to random 
excitation, using the data analysis procedures described in this report. 
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APPENDIX 


Other Data Analysis Variations of the Time-lag Products 
Correlation/Frequency Analysis Procedure 

Hanning filter . - Figures Al, A2, and A3 show the effect of applying the 
Hanning filter (or window) to the transfer functions and phase angle spectras 
for the five peak responses using FLEX excitation. Inspection of the shape of 
the functions and resulting frequency and damping lead to the conclusion that 
the Hanning filter had a negligible effect on smoothing these particular peak 
responses. 

Using the oscillator signal instead of the FLEX force . - Figure A4, a, b, 
and c, shows the effect of using the oscillator signal instead of the FLEX 
force for all five peak responses. The advantage of using the oscillator signal 
(which was used to command the FLEX wand frequency) instead of the FLEX force 
time history in the cross -correlation and cross -spectra calculations is that the 
oscillator signal has little noise, which could help smooth the spectra, while 
the force signal does suffer from noise, such as is caused by "rattle" in the 
system. The disadvantages come from the fact that the FLEX wand response 
slightly lags the oscillator command due to a finite reaction time of the 
hydraulic system and time lags introduced by rectification, which also changes 
with frequency (figure 3) . This phase lag would introduce an error in the 
cross-correlation calculations and the resulting spectra displays. Another 
disadvantage is that the oscillator signal amplitude slightly decreases with 
increasing frequency, while the FLEX force does not follow this trend, as 
shown in figure 10. This would also introduce errors in the calculations. 

Figure A4, a, b, and c, shows the gain and phase spectra using the oscil- 
lator signal (and Hanning filter). Comparison of figure A4, a, which uses the 
oscillator, with figure 22, which uses the FLEX force, shows no significant 
differences for the wing response. However, the comparison for the horizontal 
tail response in figures A4, b, and 23 shows the oscillator reduces the amount 
of definition of the gain peak, as evidenced by the longer extrapolation to the 
half-power point required by the oscillator results. This caused a slight 
increase in damping from the gain and a slight decrease in damping from the 
phase. The vertical tail/empennage response comparison in figures A4, 6, and 
24 exhibits the same characteristics as the horizontal tail. 

Thus, the use of the oscillator instead of the FLEX force slightly reduces 
the accuracy of the damping results of the timelag products method, but not to 
the extent that would invalidate the result for these cases. 
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Curve fitting, time, and frequency windowing. - Figures AS through A7 
show the least squares curve fit routine applied to the time lag products method 
for various time and frequency windows. Various order polynomials were applied 
to each mode, with the best fit of the gain and phase being shown. The wing 
mode, shown in figure A5, using wide frequency and time windows had a curve 
fit which greatly underestimated the damping but matched the peak frequency. 
Figure A5, b, shows the same mode using a 33-percent narrower frequency window 
and a 50-percent narrower time window, resulting in a fairly good curve fit. 

Figure A6, a and b, shows the curve fits for the horizontal tail mode. 

The curve fit estimated the spectra to include the entire width of all three 
peaks rather than to isolate the center peak. The windowed mode of figure A6,b 
(54-percent narrower frequency and 50 -percent narrower time window) had a 
curve fit which tended to flatten out the peak even more and gave a greater 
over-estimation of the damping. The manual frequency and damping calculation 
for the nonwindowed and windowed conditions gave essentially the same values, 
but in this instance both cases had the same frequency resolution. 

Figure A7, a, shows the vertical tail/empennage responses. The curve fit 
did not fit the spectrum due to more than one mode being present. It still 
fails when the mode is windowed by 50 -percent narrower frequency and time 
windows, as in figure A 7, b. The manual calculation shows that the effect of 
windowing has only a small effect on the frequency and damping results. Thus, 
this least squares curve fit is not presently developed to the point where it 
is a reliable tool to accurately estimate modal frequency and damping. The 
effect of reducing the analysis time windows by 50 percent and the frequency 
windows by 33 to 50 percent did not appreciably affect the frequency and 
damping results for the manual calculations . 

The data analysis results from these analysis options appear in tables I 
through V. 
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TABLE I . - WING SYMMETRIC SECOND BENDING MODE FREQUENCY AND DAMPING 

(Mach = 0.9 5» Altitude = 1,79*1 meters, q = 51*5 kN/rri 
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All frequencies in Hz. X = questionable value, not used in average, 

o = mean deviation. Sample rate = 360 samples/second. 




TABLE III. - VERTICAL TAIL/EMPENNAGE FIRST MODE FREQUENCY AND DAMPING 
(Mach = O. 95 , Altitude = 1,794 meters, q. = 51. 5 kN/m 2 ) 
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NCHHS: All frequencies in Hz X — questionable value, not used in average 

a = mean deviation sample rate = 360 samples/second 
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All frequencies in Hz X - questionable value, not used in average 

o = mean deviation Sample rate = 360 samples /second 




TABLE V. - VERTICAL TAIL/EMPENNAGE THIRD MODE FREQUENCY AND DAMPING 

(Mach = 0.95, Altitude = 1 , 79 ^ meters, q = 51*5 kN/m 

Time-lag-products autocorrelation method 
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NOTHS: All frequencies in Hz X - questionable value, not used in average 

a - mean deviation sample rate = 360 samples/sccond 













Figure 1. B-l Prototype test aircraft. 
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Figure 4a. - Flutter excitation system. 
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Figure 4b. FLEX exciter on test stand. 
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Figure 5. - Flutter excitation system control loop. 










(a) Shaker force output of FLEX wand 
with distributed mass 
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(b) Equivalent point mass 


Figure 9. - Required chordwise location of the shaker accelerometers 
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Figure 10. - Typical flutter excitation system forces input 
to structure at surface tips. 
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Figure 11. - Flight conditions used for data analysis along lines of constant 
Mach number, constant altitude, and constant dynamic pressure. 


A 1 1 i tude ( feet) 


L. 

<D 



3 

Q. 

E 

O 

u 

0) 

C7) 

<o 

u 

O 

4-1 

in 

<U 

CD 


C 

o 


CO 

c 

CO 


CO 

•o 


O) 

</* 

£ 


&0 

o 

r— I 

o 

r o 

0 

.c 

+J 

1 


to 

X 

r*H 

CO 

C 

CO 


CO 

X> 


CO 

o 


f- 

<D 


3 

4-1 

•H 

r-H 

Uh 


CD 

f-i 

3 

DO 

•<—( 

tL, 


46 






f unct i on 


a> 


0 * peak freq 


2t r 

g =* W 1 


(Dr 



(c) Co-and quadrature 
spectra 


Co .* freq = 03 0 
27T 



Quad : g 



Figure 13. - Five methods of manually calculating modal frequency 
and/or damping from the time-lag-products correlation/ 
frequency analysis procedure. 
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Figure 13- - Continued. 
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Figure 14. - General cross -spectral analysis model 
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(a) FLEX excited transfer functions (GAIN) 

using Time- 1 ag-products cross-correlation 
method 


Figure 17. - Examples of analysis results of five peak responses 
at Mach = 0.95, q = 51.5 kN/m . 
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(b) Random excitation power spectral densities of responses using 
t i me- 1 ag-product s autocore 1 1 at i on method 
Figure 17 . “ Continued. 
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Figure 18. - Time histories of five peak responses 2 
at 360 samples per second at Mach = 0.95, q = 51.5 kN/m 





Q| 

ii 


1 


II 

III 







1 

1 

II 

1 

1 





1 



i 


ii 

II 


i 

i 

i 



1 


i 

i 


1 




1 

i 



in 

Mj 

ip 





III 







R 

1 









1 



i 


K 

» 


n 

i 


i 

i 

1 


i 



1 





i 



ui 

H; 

!lll 


1 


V. 

1! 


1 





11 

1 


l 

1 

d 





1 


i 

»■ 


II 

III 


I 

n 

a 

u 

E 

a 


i 

i 


0 




1 

i 



m 

m 

I'll 


i 


VI 

I'll 







'1 

H 




1 





a 


s 

i 


31 

in 


i 

« 

k 

i 

ii 

E 


a 

£ 


ft 




a 

ft 



m 


II 




II 

III 







1 

II 



1 








i 



II 

u 



i 

■ 

i 


i 


i 

1 


1 




i 

i 



in 

pj 

II 


1 


II 

III 







1 

1 



1 

1 





i 


i 




ii 


i 


a 

i 

1 

i 


i 

1 







i 



UI 


It 17 IB l« » II 


TIME (SECONDS) 

Random response 



FLEX forced response 


lllll 

llll 

Ill 

II 

m 

hi 

n 

II 

in 

i 

II 

1 

hi 

II 

1 

llll 

in 

III 

II 

ii 

n 

II 

III 

II 

III 

r 

II1EM 

mm 

Dill 

in 

m 

HI 

m 

rn 

MI 

i 

11 

n 

in 

1 

H 

inn 

in 

im 

Ml 

in 

UI 

HU 

mil 

nn 

III 

t 

HUM 

Ml 

m 

m 

in 


Hi 

m 

mm 

Id 

31 


18 

1 

m 

ni 

881 

a 

im 

81 

IS 

m 

S3 

III 

T 

man 

rrtip 

IK 

IV 

n 



i 


R 

% 

n 

Hf 

91 

8 

mi 

m 

Bill 


« 


III 

m 

83 

III 

r 

iuiii 

llftl 

III 

ii 

ii 

i 


n 

hi 

i 

A 

ii 

HI 

II 

II 

IUI 

INI 

Lit 

UI 

ii 

Hi 

II 

IHH 

UMI 

III 

r 

- 4 — 

iiiii 



n 

n 

i 


ii 

m 


1 

ii 

II 

II 

1 

Hi 

III 

III 

ii 

ii 


II 

llll 


IU 

L 


TIME (SECONDS) 


Flutter excitation system force 



(b) Horizontal tail symmetric second bending mode 


Figure 18. - Continued. 
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(a) Autocorrelation functions 


of wing second bending response to random excitation, 
q = 51.5 kN/m^, for 25-second time history. 
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(b) Power spectra of au tocor re 1 a t i on of responses and phase 

angle of cross-cor re 1 a t i on for 36 to 66 rad/sec, frequency 
resolution = 0.25 rad/sec. 

Figure 22. - Continued. 








Figure 23. - Analysis of horizontal tail second bending 
random excitation at Mach =0.95, q = 51.5 kN/m^, for 
time history with Hanning filter. 
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(b) Power spectra of au tocor re 1 a t i ons of 
responses and phase angle of cross- 
correlation for 93 to 130 rad/sec, 
frequency resolution = 0.50 rad/sec. 


Figure 23. - Continued. 
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Vertical tail tip 



(a) Autocorrelation function 



(b) Power spectra of au tocor re 1 at i on of 
responses for 60 to 85 rad/sec, 
frequency resolution = 0.50 rad/sec. 


Figure 24. - Analysis of vertical tail/empennage response to 
random excitation at Mach =0.95, q = 51.5 kN/m 2 , 
for 25-second time history with Hanning filter. 
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(c) Power spectra of autocorrelation of 
responses for 70 to 86 rad/sec, 
frequency resolution = 0.50 rad/sec. 


Figure 24. - Concluded. 










Figure A2. - Effect of Hanning filter on horizontal tail symmetric second bending mode gain 

and phase spectra at Mach = 0.95, q = 51.5 kN/m. 







Figure A4. - Effect of using oscillator signal instead of FLEX force^on gain and 
phase (with Hanning filter) at Mach = 0.95, q = 51.5 kN/nf'. 




<t> (DEGREES) GAIN 




(c) Vertical ta i 1 /empennage modes 


Compare with figure 24, 
which uses FLEX force. 
Gain: freq = 10.70 Hz, 

9 = 

0.69, 

freq = 11.70 Hz, 

9 = 

0.0925 

freq = 13*09 Hz, 

g = 

0.073 


Phase: g = 0.0905, g = 0.147, g = ? 


Figure A4. - Concluded. 
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frequency window = 60 to 85 rad/sec frequency window = 58 to 70 rad/sec 

time window = 25 seconds time window = 12.5 seconds 

Figure A7. - Least squares curve fit of gain and phase for vertical tail/empennage modes at 
Mach = 0.95, q = 51.5 kN/m^, with various time and frequency windows with Hanning. 
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